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The  U.  S.  Air  Force,  through  a  contractual  effort  by  Argonne  National  Laboratory, 
Is  developing  a  computerized  Air  Quality  Assessment  Model  (AQAM)  to- better 
analyze  the  environmental  impact  of  airbase  operations.  While  most  emphasis 
1$  placed  on  improving  the  accuracy  of  treating  aircraft  operations,  nearly  any 
municipal  or  industrial  emission  source  can  be  analyzed.  Emission  calculations 
ere  performed  within  the  model  so  that  only  detailed  source  activity  information 
must  be  Input.  Ambient  air  quality  can  be  predicted  for  hourly,  daily,  monthly, 
or  annual  averaging  times  of  up  to  ten  pollutants.  This  report  describes  the  gen¬ 
eral  capabilities  and  structure  of  the  AQAM. 
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Development  of  an  Air  Force  Air  Quality 
Assessment  Model  (AQAM) 


D.  F.  NAUGLE 


INTRODUCTION 

The  primary  objective  of  the  U.  E.  Air 
Force  Air  Quality  Assessment  Model  (AQAM)  is  to 
predict  ground  level  ambient  air  quality  pollutant 
concentrations  which  result  from  emission  sources 
both  within  and  surrounding  Air  Force  bases. 
Special  emphasis  has  been  given  to  the  accurate 
description  and  dispersion  of  aircraft  emissions. 
However,  nearly  any  urban  point,  line,  or  area 
source  can  be  treated.  This  model  can  be  applied 
to  at  least  four  different  functions! 

1  Environmental  impact  assessments 

2  Comparison  of  predicted  concentrations 
to  state  and  federal  air  quality 
standards 

3  Prioritizing  the  significance  of  emis¬ 
sion  sources 

4  Evaluating  potential  pollution  design 
criteria  for  the  engineering  of  control 
devices . 


The  AQAM  is  being  developed  for  the  Air 
Force  by  Argonne  National  laboratory  (l).1  It 
incorporates  many  features  from  previous  urban 
dispersion  models.  Annual  pollution  levels  can 
be  predicted  using  stability  wind  frequency  dis¬ 
tributions  aB  applied  in  the  widely  used  Air 
Quality  Display  Model  (2).  Short-term  dispersion 
(1-  to  24-hr  time  averaging)  uses  the  double 
CauBsian  equations  hb  presented  in  a  dispersion 
workbook  by  D,  Bruce  Turner  (J5).  Average  urban 
emission  densities  can  be  either  directly  input 
or  calculated  from  the  average  categories  as 
used  in  the  Northern  Research  and  Engineering 
Corporation  Airport  Model  (4).  Automotive  emis¬ 
sions  are  calculated  by  using  segments  of  the 
Argonne  National  laboratory  Transportation  Model 
(jj).  All  line  sources,  including  aircraft,  use 

1  Underlined  numbers  in  parenthesis  desig¬ 
nate  References  at  end  of  paper. 
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Table  1  Emission  Sources  for  Air  Quality  Analysis 


AIRCRAFT 
t.  FIXED  WING 


GROUND  MU II 1 1  L' 

I.  AGE 


FACILITIES 

I.  POWER  PLANTS 


2.  HELICOPTER 


A.  GENERATORS  (2  TYPES) 

B.  HEATERS 

C.  COOLERS 

D.  ENGINE  STARTERS 

E.  HYDRAULIC  TEST  STANDS 

F.  LIGHT  CARDS 

G.  COMPRESSORS 

H.  TURBINE  AGE  EQUIPMENT 
|.  OTHER 

2.  LIGHT  DUTY  VEHICLES 

MILITARY 

3.  LIGHT  DUTY  VEHICLES 
CIVILIAN 

4.  TRUCKS 

5.  OFF  ROAD  VEHICLES 


2.  INCINERATORS  (PATHOLOGICAL) 

3.  INCINERATORS  (WOOD  &  PAPER) 

4.  AIRCRAFT  FUELING  FACILITIES, 

5.  CAR  l  TRUCK  FUELING  FACILITY 

6.  TANK  TRUCK  LOADING  UNLOADING. 

7.  DISTRIBUTION  TANK 

8.  STORAGE  TANKS 

9.  HEATING  PLANTS 

A.  AIRMEN  DORMITORY 

B.  BOO 

C.  FAMILY  HOUSING 

10.  TEST  CELLS 

11.  RUN  UP  STANDS 

12.  INDUSTRIAL  FACILITIES 

13.  TRAINING  FIRES 


the  Argonne  Line  Source  Model  (6), 

Unique  features  of  the  A  QAM  Include  very 
detailed  algorithms  to  accurately  describe  air¬ 
craft  operations,  Aircraft  emissions  are  not 
grouped  Into  generalized  categories  (e.g,,  air¬ 
craft  with  over  400,000  lb  gross  weight).  In¬ 
stead,  engine  mission  factors,  aircraft  landing 
and  take-off  (LTO T  cycle  times,  runway  roll 
distances,  take-off  profiles,  and  approach  pro¬ 
files  are  all  treated  specifically  for  the 
particular  aircraft  types  of  interest.  Also, 
atmospheric  conditions  are  treated  in  more  de¬ 
tail  by  the  AQAM  than  by  previous  long-term 
dispersion  models. 

GENERAL  DESCRIPTION  OP  MODEL 

The  A QAM  presently  has  four  functional 
elements  as  presented  In  Fig,  1.  The  Source 
Inventory  Program  accepts  operational  input  data 
and  computes  the  true  annual  mean  emissions  for 
each  source.  This  output  alone  is  very  useful 
and  may  be  adequate  for  many  environmental  Impact 
assessments  without  any  further  dispersion  analy¬ 
sis,  When  dispersion  calculations  are  needed, 
tills  true  annual  mean  emission  is  adjusted  to  the 


actual  monthly,  weekly,  and  dally  time  period 
emission  rate  for  use  in  the  Short-Term  or  Long- 
Term  Dispersion  Model, 

Operational  input  to  the  Source  Inventory 
Program  1b  broken  into  three  categories  as  shown 
in  Pig.  1.  AIRCRAFT  DATA  includes  all  emissions 
directly  related  to  operations  of  the  aircraft. 
Aircraft  taxi,  landings,  take-offs,  fuel  venting, 
and  ground  support  equipment  are  examples,  AIR¬ 
BASE  DATA  includes  all  emission  sources  within 
the  confines  of  the  airbase  which  are  in  the 
aircraft  operations  category.  Examples  are  Jet 
engine  test  cells,  power  plants,  heating  plants, 
military  industrial  operations,  vehicular  opera¬ 
tions,  and  space  heating  of  buildings,  ENVIRONS 
are  defined  as  all  emission  sources  outside  the 
airbase.  Examples  are  municipal  power  plants, 
civilian  industry,  and  civilian  vehicles. 

The  Short-Term  Dispersion  Program  calculates 
air  quality  by  using  hour  by  hour  meteorological 
input  data.  Computations  can  be  averaged  for 
time  periods  up  to  24  hr.  This  allows  direct 
comparisons  to  federal  ambient  air  quality 
Btandarda  (which  may  be  specified  for  a  1-,  3~, 

8-,  or  24-hr  averaging  basia,  depending  on  the 
pollutant  type). 


Table  2  Source  Inventory  Output  for  Aircraft 
T  33  AIRCRAFT  -  ANIIUAL  MISSIONS  (METRIC  T0N8/YRAR) 


OPERATION 

Cl 

HC 

NOX 

PM 

SOX 

STARTUP 

r.6  ‘kit  *i. 

1.3 552*92 

2.tl3c*Cl 

1.139E4S2 

1.5  19l+U3 

TAXI  OUT 

1.  c».i  :*i  3 

Z.OibL*  32 

3. 6  1  b E  ♦  0  1 

2.161E+J2 

2.«62t>30 

ENGINE  CM:C< 

7 ;  c :  ♦  c  1 

5.927E03 

5.318L+CC 

1. 3  22E*3  1 

4.3 87t“ 31 

RUNWAY  R3LL 

e.  4  1 

9 . 0  7  7  E  ♦  G  3 

7.5i!:i>Ji 

1.21*i»31 

9.677E-IH 

CL  I'll  1 

7.0'-?i*ii 

s.:F«.F*ac 

8.2  42E.  +  V  C 

l.ij7;+vi 

8 .  k  54c-iii 

CLI'11  2 

1 . 4 «.  ■: :  *  c  2 

1 . 9 1 5 :  ♦  G 1 

1.69kfc01 

3.l3JE»3i 

1. 1 52 £♦ 0  & 

APPROACH  i 

1  •  »  b£  ♦  C  t 

2.uC3£»Cl 

9.97C8+QC 

5. 5  3  9£ ♦  J  1 

1.138tti)i) 

APPROACH  ? 

1.  2-2CU  2 

2 . 0  7  9£ ♦ t  1 

£..  692tOC 

3«7Jb£t"ii 

b.  8 llc-u  1 

LA  NO INC 

1.3.  '<.>.2 

3  .  rt  C  :  t  C  1 

4.831E*9C 

*.71k£*]1 

3.b24t-'Jl 

TAXI  IN 

A.  6t'4i»C2 

2.*33-:*32 

3.97CL»‘31 

1. 7 3BE »  1 2 

2.  J l7_c  l  J 

SHUTDOWN 

?. 3-o:*Cl 

1 .  975£  ♦ C  0 

1. OCoL+QC 

5. 5  9 0 i  ♦  3  C 

7.59vl-32 

ARP  ♦  QE»  SV 

7.  2  tj 

1.15JEK  2 

2.  <» 7 3 £ ♦  C  L 

1.191  20  1 

4.  b  J7l*1  1 

FUEL  VfNTING 

*  < 

1 . 195C ♦ 0  2 

C. 

C. 

C. 

FILL  ♦  VIU 

.  • 

2.551.E*  )2 

a. 

C. 

c. 

T3UC1  ♦  G1 

7.  5c  5; »  3? 

1. 1  d  3  £  ♦  0  J 

5. 976E»3 1 

1.9S2E»02 

5.8J9c*QQ 

TOTAL 

4.  3  n  3£  »  C  3 

1  «»73E* 13 

2,G99£»Q2 

9.1 J2i*  3 2 

1.888L»U1 

TRANSENT 


Table  3  Souroe  Inventory  Summary  Output 
SUMMARY  OF  ANNUAL  EMISSIONS  (METRIC  TONS/YEAR) 


0PGRATI1N 

Cl 

HC 

NOX 

PM 

SOX 

A  I °CRAFT 
GROUND  M13TL 
FACILITIES 
ENVIRONS 

7.9  j  &  £  ♦  C  3 
?,  5r7£ ♦  03 

5.  629EKI 
*\ 

*• 

2.6C3E  +  C  3 
J.997E02 
l.-65E*C2 

j  « 

4.39bc+)2 

3.  49bL  +  32 

3. 34i£* J  2 

C. 

1.7  811*3 3 
2.7 1  i  £♦ C  1 

1.7 ]«£♦ 12 

C. 

3.694t»dl 
1. UbottOl 

3*247c»00 

C  • 

GRAND  TOTAL 

1. 125EK4 

3. 149E+C  3 

1.393E»03 

1.179E03 

5.383£»0l 

PERCENT  OF  EMISSIONS  FROM  ALL 

SOURCES 

opiratiin 

Cl 

HC 

NOX 

PM 

SOX 

RIPCRAFT 
GROUNO  M01TL 
FACILITIES 
ENVIRONS 

C3.397 

31.898 

5.3(5 

O.OCO 

82.  ‘>5 4 

12. 393 
4.  *53 
C.  303 

43.238 

31.97fc 

27.918 

0.300 

91.301 

1.370 

d.oJu 

0.030 

88 » 0  3  7 
25.330 
b.C  33 

o.ooc 

The  Long-Term  Dispersion  Program  utilizes  consists  of  specialized  climatological  data 

extensive  climatological  data  from  the  Meteoro-  categorized  by  diurnal  time  periods  on  a  monthly 

logical  Data  Program  and  emissions  from  the  and  annual  basis.  This  is  a  unique  feature  of 

Source  Inventory  Program  to  produce  monthly,  AQAM  which  is  presently  undergoing  a  value 

seasonal,  or  annual  average  ambient  concentre-  analysis, 

tlons.  Averages  can  also  be  performed  for  one 

of  several  daily  time  periodB  (e.g,,  annual  SOURCE  INVENTORY  PROORAM 

average  air  quality  for  the  time  period  0600  to 

0900  can  be  calculated).  Analysis  can  be  done  Manually  calculating  a  total  emission  in- 

for  all  pollutants,  but  is  especially  important  ventory  can  be  a  tedious  procedure.  This  is 

for  oxides  of  Nitrogen  which  presently  only  have  particularly  true  for  aircraft  which  have  highly 

an  annual  ambient  air  quality  standard.  variable  engine  emission  characteristics  and  many 

The  Meteorological  Data  Program  produces  engine  modes  throughout  the  landing  and  take-off 

e  meteorological  data  tape  (for  input  to  the  cycle.  This  Source  Inventory  Program  was.  there- 

Icng-Term  Dispersion  Program).  The  data  tape  fore,  designed  to  contain  data  banks  of  emission 
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information.  Total  emissions  are  then  calculated 
when  given  only  source  activity  information  from 
a  field  survey.  Jet  engine  test  cell  emissions, 
for  example,  can  be  calculated  by  knowing  only 
the  type  of  aircraft  engines  tested  and  the 
operating  times  in  each  engine  mode. 

The  Holland  (J)  and  Carson-Moses  (8)  plume 
rise  relationships  are  programmed  into  the  com¬ 
puter  code  for  optional  usage.  Either  can  be 
used  with  input  information  of  stack  height, 

Btack  exit  temperature,  gas  exit  velocity,  stack 
diameter,  and  building  height. 

Examples  of  typical  emission  sources  on  an 
Air  force  Base  are  shown  in  Table  1.  Aircraft 
are  given  an  area  source  geometry  at  the  parking 
ramp,  and  line  geometries  for  taxi,  take-off, 
climb-out,  approach,  and  landing.  Coordinates 
of  all  area  and  line  sources  must  be  input. 
Aerospace  Oround  Equipment  (AQE)  are  considered 
area  sources  while  vehicles  can  have  either  a 
line  or  area  source  geometry,  facilities  can  be 
treated  as  having  point,  line,  or  area  Bource 
emission  geometries.  Also,  any  sources  not 
specifically  programmed  can  be  added  without 
computer  code  changes  by  the  use  of  NAMEUST 
in  pit  statements , 

Sample  source  inventory  computer  outputs 
for  aircraft  is  presented  in  Table  2,  Similar 
outputs  are  generated  for  each  aircraft  type  and 
for  totals  of  all  aircraft.  All  components  re¬ 
lated  to  the  aircraft  LTO  cycle  are  listed 
separately  so  that  the  effects  of  any  pollution 
control  strategies  can  be  readily  analyzed. 

Similar  types  of  summary  charts  are  also  auto¬ 
matically  printed  for  all  AIRBASE  and  ENVIRON 
emission  sources. 

Computer  output  summary  tables  of  all  pre¬ 
vious  source  calculations  are  shown  In  Table  3. 
Total  annual^emisalons  are  given  for  up  to  ten 
poll\jfc&tttiTof  Interest,  The  percent  of  each 
.etfurce  category  toward  the  total  Is  bIbo  calcu¬ 
lated. 

SHOHT-TERM  DISPERSION  PROGRAM 

The  Short-Term  Dispersion  Program  is  used 
to  predict  air  quality  concentrations  resulting 
from  specifically  input  neteorological  conditions. 
These  predictions  are  important  when  assessing 
the  Impact  of  sources  releasing  acutely  toxic 
substances,  such  as  carbon  monoxide,  or  when 
comparing  predicted  results  with  limited  ambient 
air  quality  measurements. 

The  first  task  of  this  program  is  to  read 
in  the  Source  Inventory  tape  containing  the  cal¬ 
culated  true  annual  mean  emissions  for  each 
souroe.  These  annual  mean  emissions  are  then 
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Pig.  2  Example  of  activity  pattern  Information 


adjusted  to  the  specific  time  period  of  interest. 
For  example,  most  military  aircraft  have  decreased 
flying  activities  during  weekends  and  after  dark¬ 
ness.  Inaccuracies  would  result  from  simply 
dividing  the  true  annual  emissions  by  the  total 
number  of  hourB  per  year  to  get  an  hourly  emis¬ 
sion  rate.  Activity  patterns  as  graphically 
shown  in  Pig,  2  are,  therefore,  input  into  the 
Short-Term  Program  and  emissions  are  automatically 
calculated  for  the  particular  month,  day,  and 
hour  of  interest.  By  programming  aotlvity  pat¬ 
terns,  which  are  usually  available  from  paat 
records,  the  logistical  problem  of  recalculating 
total  emissions  for  each  time  period  is  avoided. 
Without  activity  patterns,  emissions  rates  would 
have  to  be  broadly  categorized  (e.g.,  emissions 
per  weekday).  This  would  minimize  calculations 
but  bIbo  Introduce  inaccuracies. 

The  next  task  of  the  Shogt-Term  Program  Is 
to  process  the  hourly  Input  of  specific  meteoro¬ 
logical  conditions.  Atmospheric  stabilities  (.2), 
wind  speed,  wind  direction,  temperature,  and 
mixing  depth  are  input  and  can  be  changed  on  an 
hourly  basis. 
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Dispersion  calculations  can  now  ba  done 
since  the  speolfic  source  emission  and  dispersion 
conditions  are  known.  Computations  are  done  for 
each  of  several  thousand  possible  sources,  and 
the  results  are  summed  at  each  of  30P  possible 
receptors.  Ambient  concentrations  at  each  re* 
ceptor  are  calculated  for  three  categories  — 
AIRCRAFT,  AIRBASE,  and  ENVIRONS  —  so  that  pol¬ 
lution  contributions  tor.  each  category  can  be 
identified. 

l^e  final  task  of  the  Short-Term  Program 
is  to  output  the  data  in  a  usable  format.  Re¬ 
sults  can  be  printed  hourly  or  averaged  over 
3-,  8-,  or  24-hr  periods,  'a  example  o i  the 
tabular  output  ip  shown  ih  Table  4.  Similar 


tables  are  produced  for  AIRCRAFT,  AIRBASE, 
ENVIRONS,  and  TOTAL,  Results  can  also  be  plotted 
for  graphical  presentation, 

LONG-TERM  DISPERSION  PROGRAM 

The  Long-Term  Program  performs  many  more 
calculations  than  the  Short-Term  Program.  Cal¬ 
culations  to  determine  the  contribution  of  each 
source  to  the  many  receptors  must  be  done  sep¬ 
arately  for  16  wind  directions,  Bix  wind  speed 
classes,  and  five  atmospheric  stability  condi¬ 
tions.  Well  over  1  billion  computations  can  be 
required.  These  are  all  in  addition  to  the  data 
storage,  data  conversion,  summation  routines, 
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and  other  "bookkeeping”  procedures  to  produce  the 
output  in  a  meaningful  format,  A  centra:  Memory 
Unit  Core  Storage  capacity  of  300,000.  60-bit 
worde  on  the  CDC  6600  computer  are  typically 
used. 

Several  user  options  can  be  utilized  in 
the  Long-Term  Program  so  that  adequate  analyses 
can  be  performed  but  unnecessary  computations 
are  avoided.  The  simplest  case  is  to  request 
only  the  arithmetic  annual  mean  concentrations 
for  the  total  24-hr  period  within  each  day.  If 
more  specific  predictions  are  desired,  they  can 
be  made  separately  for  each  month  of  the  year. 
Diurnal  patterns  can  also  be  readily  analyzed  by 
making  monthly  predictions,  in  any  of  the  follow¬ 
ing  dally  time  periods i 

1  0000  thru  2400  (all  day) 

2  0600  thru  1800  (daylight  hours) 

3  0600  thru  0900  (peak  traffic  hours) 

4  0900  thru  I5OO  (major  working  hours) 

5  1500  thru  1800  (afternoon  traffic  hours) 

6  1800  thru  2100  (nightly  inversion  forma¬ 
tion) 

7  2100  thru  0C00  (overnight  hours). 

Only  the  specific  meteorological  conditions  and 
emission  rates  will  be  used  within  each  of  the 
foregoing  time  periods  of  interest. 

The  ability  of  the  Long-Term  Dispersion 
Program  to  utilize  the  diurnal  patterns ‘of  both 
meterology  and  emissions  is  another  unique  feature 
of  A QAM,  An  annual  or  monthly  average  can  be 
made  using  time  periods  3  thru  7.  An  average 
formed  in  this  maimer  should  give  results  con¬ 
siderably  closer  to  reality  than  one  using  daily 
average  emissions  and  meteorological  conditions. 
For  example  consider  the  fact  that  ambient  con¬ 
centrations  (X)  result  in  the  atmosphere  by  the 
dispersion  of  emissions  (Q).  Therefore, 

X  -  f (Q)  (1) 

The  atmosphere's  dispersive  ability  la  represented 
by  the  function  "f."  Since  both  "Q"  and  ”f”  are 
strongly  dependent  on  the  time  of  day,  the  re¬ 
sulting  mean  concentration  (x)  for  a  daily  period 
is  more  accurately  calculated  byi 

^per  day  "  ?  fW  4  F2(Q2>  *  —  Wil  <2) 

a 

Rather  than  the  commonly  used  method  oft 

f  7  (Q) 

Aper  day  ■  *per  day  per  day  (3) 


The  accuracy  Improvements  by  using  this 
technique  have  not  yet  b?en  quantified.  Argonne 
National  laboratory  is  presently  Investigating 
this  difference  using  A QAM,  They  are  comparing 
long-term  air  quality  predictions  in  the  Chicago 
area  using  the  five  dally  time  periods  to  those 
made  using  single  daily  means  of  emissions  and 
meteorological  conditions.  The  results  will  than 
be  compared  to  aotual  air  quality  data. 

The  output  of  the  Long-Term  Dispersion 
Program  is  in  exactly  the  6ame  format  as  for  the 
Short-Term  Program,  Isopleth  contours  (equal 
pollution  lines)  of  the  Long-Term  results  art 
especially  significant  because  winds  from  all 
directions  have  been  considered.  The  resulting 
isopleths  often  have  considerably  different 
patterns  than  those  of  a  wind  rose.  This  is  due 
to  the  effects  of  different  atmospheric  stabili¬ 
ties  and  wind  speeds  on  the  dispersion  of  pollu¬ 
tion  sources  whioh  are  elevated  above  ground 
level. 
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METEOROIOOICAL  DATA  PRO  OR  AM 

To  supply  climatological  information  re¬ 
quired  by  the  Long-Term  Dispersion  Program,  a 
computer  tape  of  meteorological  data  must  be 
compiled  specifically  for  eaoh  location  of 
Interest,  These  tapes  are  compiled  from  data 
which  is  stored  at  the  National  Climatic  Center, 
National  Oceanic  and  Atmospheric  Administration 
(NOAA),  Asheville,  North  Carolina,  Past  records 
of  periods  between  5  and  20  years  are  typically 
available.  The  Meteorological  Data  Program 
described  in  this  section  is  being  operated  for 
Air  Force  agencies  by  the  USAF  Environmental 
Teclmical  Applications  Center  (ETAC),  Washington, 
DC. 

Information  la  structured  to  first  give  the 
average  annual  conditions  as  follows: 

1  Average  temperature  (deg  F) 

2  Average  number  of  degree  days  (days  • 
deg  F) 

3  Average  pressure  (millibars) 

4  Average  pressure  altitude  (100  ft) 

5  Average  wind  speed  (m/seo) 

6  Average  daily  temperature  variation 
(deg  F) 

7  Average  afternoon  maximum  depth  (m). 

Data  is  then  categorized  into  daily  time 
intervals  of:  0000  thru  2400,  0600  thru  1800, 
0600  thru  0900,  0900  thru  I500,  1500  thru  1800, 
1800  thru  2100,  2100  thru  0600.  Within  each  of 
these  time  periods,  data  is  stratified  according 
to  the  following  scheme  in  Fig.  3.  The  widely 
used  Martin-Tikvart  frequency  distributions  (10) 
compute  average  mixing  depths  over  the  entire 
period.  Greater  accuracy  should  be  obtained 
by  using  this  array  which  varies  mixing  depth  as 
a  function  of  wind  speed,  stability,  diurnal 
time  period,  and  month.  The  degree  of  this  ac¬ 
curacy  improvement  la  currently  being  studied  in 
a  sensitivity  analysis  of  the  model. 
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